142 research outputs found
On a Cahn--Hilliard--Darcy system for tumour growth with solution dependent source terms
We study the existence of weak solutions to a mixture model for tumour growth
that consists of a Cahn--Hilliard--Darcy system coupled with an elliptic
reaction-diffusion equation. The Darcy law gives rise to an elliptic equation
for the pressure that is coupled to the convective Cahn--Hilliard equation
through convective and source terms. Both Dirichlet and Robin boundary
conditions are considered for the pressure variable, which allows for the
source terms to be dependent on the solution variables.Comment: 18 pages, changed proof from fixed point argument to Galerkin
approximatio
Long time existence of solutions to an elastic flow of networks
The L2-gradient flow of the elastic energy of networks leads to a Willmore type evolution law with non-trivial nonlinear boundary conditions. We show local in time existence and uniqueness for this elastic flow of networks in a Sobolev space setting under natural boundary conditions. In addition, we show a regularisation property and geometric existence and uniqueness. The main result is a long time existence result using energy methods
Efficient cosmological parameter sampling using sparse grids
We present a novel method to significantly speed up cosmological parameter
sampling. The method relies on constructing an interpolation of the
CMB-log-likelihood based on sparse grids, which is used as a shortcut for the
likelihood-evaluation. We obtain excellent results over a large region in
parameter space, comprising about 25 log-likelihoods around the peak, and we
reproduce the one-dimensional projections of the likelihood almost perfectly.
In speed and accuracy, our technique is competitive to existing approaches to
accelerate parameter estimation based on polynomial interpolation or neural
networks, while having some advantages over them. In our method, there is no
danger of creating unphysical wiggles as it can be the case for polynomial fits
of a high degree. Furthermore, we do not require a long training time as for
neural networks, but the construction of the interpolation is determined by the
time it takes to evaluate the likelihood at the sampling points, which can be
parallelised to an arbitrary degree. Our approach is completely general, and it
can adaptively exploit the properties of the underlying function. We can thus
apply it to any problem where an accurate interpolation of a function is
needed.Comment: Submitted to MNRAS, 13 pages, 13 figure
A Phase Field Model for Continuous Clustering on Vector Fields
A new method for the simplification of flow fields is presented. It is based on continuous clustering. A well-known physical clustering model, the Cahn Hilliard model, which describes phase separation, is modified to reflect the properties of the data to be visualized. Clusters are defined implicitly as connected components of the positivity set of a density function. An evolution equation for this function is obtained as a suitable gradient flow of an underlying anisotropic energy functional. Here, time serves as the scale parameter. The evolution is characterized by a successive coarsening of patterns-the actual clustering-during which the underlying simulation data specifies preferable pattern boundaries. We introduce specific physical quantities in the simulation to control the shape, orientation and distribution of the clusters as a function of the underlying flow field. In addition, the model is expanded, involving elastic effects. In the early stages of the evolution shear layer type representation of the flow field can thereby be generated, whereas, for later stages, the distribution of clusters can be influenced. Furthermore, we incorporate upwind ideas to give the clusters an oriented drop-shaped appearance. Here, we discuss the applicability of this new type of approach mainly for flow fields, where the cluster energy penalizes cross streamline boundaries. However, the method also carries provisions for other fields as well. The clusters can be displayed directly as a flow texture. Alternatively, the clusters can be visualized by iconic representations, which are positioned by using a skeletonization algorithm.
Mean curvature flow with triple junctions in higher space dimensions
We consider mean curvature flow of n-dimensional surface clusters. At
(n-1)-dimensional triple junctions an angle condition is required which in the
symmetric case reduces to the well-known 120 degree angle condition. Using a
novel parametrization of evolving surface clusters and a new existence and
regularity approach for parabolic equations on surface clusters we show local
well-posedness by a contraction argument in parabolic Hoelder spaces.Comment: 31 pages, 2 figure
Smolyak's algorithm: A powerful black box for the acceleration of scientific computations
We provide a general discussion of Smolyak's algorithm for the acceleration
of scientific computations. The algorithm first appeared in Smolyak's work on
multidimensional integration and interpolation. Since then, it has been
generalized in multiple directions and has been associated with the keywords:
sparse grids, hyperbolic cross approximation, combination technique, and
multilevel methods. Variants of Smolyak's algorithm have been employed in the
computation of high-dimensional integrals in finance, chemistry, and physics,
in the numerical solution of partial and stochastic differential equations, and
in uncertainty quantification. Motivated by this broad and ever-increasing
range of applications, we describe a general framework that summarizes
fundamental results and assumptions in a concise application-independent
manner
Causal deep learning models for studying the Earth system
Earth is a complex non-linear dynamical system. Despite decades of research and considerable scientific and methodological progress, many processes and relations between Earth system variables remain poorly understood. Current approaches for studying relations in the Earth system rely either on numerical simulations or statistical approaches. However, there are several inherent limitations to existing approaches, including high computational costs, uncertainties in numerical models, strong assumptions about linearity or locality, and the fallacy of correlation and causality. Here, we propose a novel methodology combining deep learning (DL) and principles of causality research in an attempt to overcome these limitations. On the one hand, we employ the recent idea of training and analyzing DL models to gain new scientific insights into relations between input and target variables. On the other hand, we use the fact that a statistical model learns the causal effect of an input variable on a target variable if suitable additional input variables are included. As an illustrative example, we apply the methodology to study soil-moisture–precipitation coupling in ERA5 climate reanalysis data across Europe. We demonstrate that, harnessing the great power and flexibility of DL models, the proposed methodology may yield new scientific insights into complex non-linear and non-local coupling mechanisms in the Earth system.</p
On a diffuse interface model for tumour growth with non-local interactions and degenerate mobilities
We study a non-local variant of a diffuse interface model proposed by
Hawkins--Darrud et al. (2012) for tumour growth in the presence of a chemical
species acting as nutrient. The system consists of a Cahn--Hilliard equation
coupled to a reaction-diffusion equation. For non-degenerate mobilities and
smooth potentials, we derive well-posedness results, which are the non-local
analogue of those obtained in Frigeri et al. (European J. Appl. Math. 2015).
Furthermore, we establish existence of weak solutions for the case of
degenerate mobilities and singular potentials, which serves to confine the
order parameter to its physically relevant interval. Due to the non-local
nature of the equations, under additional assumptions continuous dependence on
initial data can also be shown.Comment: 28 page
Model order reduction approaches for infinite horizon optimal control problems via the HJB equation
We investigate feedback control for infinite horizon optimal control problems
for partial differential equations. The method is based on the coupling between
Hamilton-Jacobi-Bellman (HJB) equations and model reduction techniques. It is
well-known that HJB equations suffer the so called curse of dimensionality and,
therefore, a reduction of the dimension of the system is mandatory. In this
report we focus on the infinite horizon optimal control problem with quadratic
cost functionals. We compare several model reduction methods such as Proper
Orthogonal Decomposition, Balanced Truncation and a new algebraic Riccati
equation based approach. Finally, we present numerical examples and discuss
several features of the different methods analyzing advantages and
disadvantages of the reduction methods
A Cahn-Hilliard-Darcy model for tumour growth with chemotaxis and active transport
Using basic thermodynamical principles we derive a Cahn--Hilliard--Darcy model for tumour growth including nutrient diffusion, chemotaxis, active transport, adhesion, apoptosis and proliferation. The model generalise earlier models and in particular include active transport mechanisms which ensures thermodynamical consistency. We perform a formally matched asymptotic expansion and develop several sharp interface models. Some of them are classical and some new ones which for example include a jump in the nutrient density at the interface. A linear stability analysis for a growing nucleus is performed and in particular the role of the new active transport term is analysed. Numerical computations are performed to study the influence of the active transport term for specific growth scenarios
- …